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1. IBTaOOUGtlOl 


Fiew separation is a cowwon ©nginasrJjig pnasoisanon* 

It *»?ay inewas® or iaeroasa the maefninass of an ©nginaariiig 
isflea* Many flniiie ioiricas are based on flow separation 
Tidiil© the stall of an airfoil at large arigle of Incidence 
results in large drag and loss of lift* Flow separation 
can be caused by an adverse pressure gradient idiicli retards 
the fluid, the effect being ’^ore pronounced on the low 
!?iO«?entu» fluid near the wall* fhis fluid ultimately unable 
to overcome the ©ppostog pressure ieHeets into the streiw 
causing a reversal of flow inwediately this 

is termed •boundary layer separation*. Fertiaps the best 
knowa example of such behariour is the flow separation from 
the rear of a circular cylinder placed transverse to a 
stream of fluid* 

' In ^aglneering practice a common cause of flow 
separation Is an abrupt change In the profile of a solid 
surface in contact wil* the fluid* In an unconfined flow 
it may be caused by a depretsicm in, or a protrusion frum, 
an etherwlte smooth surface , such as a step in an ©pen channel* 
An example in confined flow 'is the change in diiaieter of a 
pipe Of tine Change In width of a channel* 



m 


. I&witigatioBs on flov sepaaraiioia haw pfoc««a«d 
ir* sawiaX &lffemn% diiractloKs* On# class ha* Saalt with tti# 
Question of 'Whether a toimdaj^ layer separates mOer a glwn. 
pressuwi distribution if so where is the ieparat.lon point* 
Another ©lass deal* with the in vise id analysis of a separated 
flow past a body in iidiich the separating strea’^l'ine Is 
specified by somw geosietric or dynamic condition <Fw@ 
streamline theories)* This class predicts drag a* a function 
of a base pressure coefficient* Still another class deals 
with flows in carlties* This class has the a.Svmt&ge that 
the length seales .in the two directions are . determined by 
the geometry of the solid boundaries and the behawiour at 
large Reynolds number is such that a core »gion of unlfow! 
worticity develops* Batchelor Cl@d®) first pointed out this 
behaviour and it has received support from direct numerlcail 
computations and eacperimental data* 

Although ths' lumber of Investigations of separated 
flows is extremely large t it cannot be said that a definitive 
basis has been laid down* iHoshko)* Even direct numerical 
integ«ition of Eavler^tokes oQuations encounters difficulty 
as laynolds number becomes large C e*g* Hung and Mataino^ 1966)* 
Xn the recent past there has be^ a trend to use the asye^totic 
me^d to handle the flow situations idiich otherwise looked 
to be intractable* This investigation is taken up with 
a similar intentim* 





of til® t>oimdiiiT* JUlaOt utiffisric&l rosults of Emg aafi Macagno 
Can l>e tiaoi to cbocli tb® a|»pfOp*lat®fioss of tli® 


axpansions* 

Mban separation oeenrs in pip® or ebenuiel aysta^y or 
in rotor t stator or casing of pmpBf turbine* or c<Mpr®*«or»t 
the resulting flow is a mmbor of tli« class internal ssparatoS 
flow* ^^«re ar® two waior f@atur«s of tbi* class* ©ft«a 
sorticity bas diffuso^l #wr ’eost of tb« fluii, unlUc® 
wbat bappans in «ict«mal flow* Mso continuity ani confina* 
wiant of tba flow by tha wall proiucas m intaraction batwasu 
tba ragion closa to tba wall and ration sway from tba wall 
and ibis intaraction strongly influancas prassura distribution* 
In astamal flows tba prassura distributiMi is by and largo 
datarminad by tba flaw away from tba wall* 

fba matbamatiaal problam will ba to solwa tba go^rnint 
aguatiM idiicb is alliptia in natura« witb prascribad muidltion 




Certain conditions sneli as fully developed 
flow at totti the sections are used and these taken 

to describe the actual conditions approxl»^atelj* But then 
the task is not si«’pl© as the theory of such paiPtial iifferen** 
ttal equations is in a far fro’w satisfactory state* Instead 
the us® is «ad@ of the asymptotic theory by taking a liwit 
in irtiich leysold number approaches infinity*, fhis limit 
is mainly desigaed to take fullest advantage of large Beynolds 
numfcers often found in practice* It provides a systematic 
method of constructing successively «ore refined approximations. 
Also it points out the domlnai^t forces in different refioni* 

Based on the computation performed by Wmg and Hacagno 
Cl®665, a simplification is possible* It is argued in an 
all together a different way that the govemlnf e^tlon 
Cfor the most of the flow regions) are like boundaapy layer 
equations* Itiis peifiaps supports the intutive feeling ^at 
viscous forces and inertia forces are Important In the Imiest 
order appraximation and mostly the information will be 
conveyed in the ^temstresa direction* this approach and 
ether aspects of formulation*©! the problem and outline of 
the method are discussed in Chapter S.# 

A new family of viscous layers termed’ weah mixing 

fe' , 

laysr were idmtified at early stages of the pressnt 
investigation* Veak. miMMg layers are visceits layers across ' 
which shasf changes In vovtlcity tidees place but not in 



X, the last chaptar iaals with 
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2*1 ASsyMPTiois oi‘ mi phsbbm i 

Floif thJpomgli a sudden expansioii Is eouslds'red* 

Th® flow is assu’^ed to Be steady, tw©-dl«!ecs tonal ImtnmT 
flow of an incowpwsslfcle Bewtonioia fluid of oonstaiit 
pfopertles* 

me fealf-widtlj of the olittmel after the expansion 
and th® airerage velocity after the expansion is used as 
length and velocity scales for iioii-dl»iienslonalising* Hie 
govemlRf equations then a^ 

**,-0 

♦ nty • - p, * a"* ♦ Mjy) (a.l) 

♦ Wy • • fy e !'■ * I%|| ♦ l^y) 

i^re X is the^ (divensionlessl downstreasi eo^ordinate 
and f is the Cdinensionless) transverse eo-ordinatei 
the origin heing at the expsnsion in the centre of ehannel* 
m and v are the non«*di»!^isionalised veloeities in and 
y<k4iteetiofui respeotively and p is the non^dimensiMalieed 
pressure* 1 is the heynold ni»her Based on half channel 
width and average velocity after the cxpension* the suffixes 
diNuite partial differantiati^* 


fli® ii»If -width of ehmml \mtom mpmslou 

is domotod ty a* 

Tfe« ‘booudary oooditioos are 


3C » • 0 1 

m « CS/2al [l-(y/al^J 

0<|y/a|< 1 



m 0 

o<SL^iri< 1 



V * 0 

0 < jru i 


**«?► 00 1 

{2/2} (1-y^l 

in ^ I 

<2*2tl 


©*»©* © 

0 ^ ly ( ^ 1 


y « t i : 

at « w *t 0 

© ^ X < 00 

(@»2e) 

til© flow 

la fmiftlier assumed to to sy«iw®trical afeomt 

a-aala* 




y » 0 : 

Wy » Of w » 0 

. © 4 * < 00 

|S»8i| 

daly tlio 

wppor half of tho ehanool nood to eoaaidoyoS 

in tho aoalyais 

o» aecoToat of syfwotpy* 
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Sow wo omsidor tlio fiow at iargo Itoyaold nwotow* ’ 

It ia fait ttiat tho fluom niqr to dowidod into oonoi^ 
fogima taaod oo aomo psofowtloa oharaotoxlatio of tto 
fogi<ma« flia«i oa^i logioii oan handlod aopamtoly tsal 
lator natotitog in a Itottiiig aoaao em to iioffoMii* 

Jftor tlio awddoa mpmBiim thaoo diffofoot togiona 
oiA to Mootifioiy tim m&SM flo«i th» tiittio m$ l^O'slxiiig 

■:f' ■ ■ 

l4iyo3F* fli* «ialB flow foMO th» mm and axp'Oaia to fill 




iUiother Ictoiis for tills solootlon of loogtlt soaios 
esn. «a.i») Im gtvm* liot tlu» longtli sc«3.«s Im x* and 
Slrtetions bo Hi and rospaotiiroly* fhm tho ioortla 




r«»pe«tiir«Xy* th.9 has th® se’??« laiigth scal®» . 

as th« eor®« Since the mss flov is small JJi th® mixing . 
layer, th® net mass flow in th® 'bufctl© is of the saw® ori«r 
as in th® oore* lenc® the velocity scales are also th® 
saw®* 

Sine® th® mixing layer has a velocity fi@li Iniwced 
hy the cor® and since it entrains th® fluid from th® huhhl®, 
its velocity scales are th® same as th® cor® and th® bufchl®* 
Further its length scale 3n x*dlr®ction is clearly th® ssms 
as the core* It th«n follows from th® continuity that the 
length seal® In y-4ir«cti««i also must h« th® same as the cor®# 

M® therefor® introduce th® new variables 

X m x/M t T * f 

0 • « I ■ ? * -ihr 

fh® following options are then ohtainet 
® ® 

yw ♦ R' ' C Wjr ♦ FFw ♦ Vj|w- 

m now consider th® limit m for fix§it X and f ♦ 
fh® assumption based on th® physical reasoning of the earlier 
paragraphs is that the limit m for fined X«f doscrihes 
the hehawiimr In all the three regions* the llmttlaif expiations 


(@•33 


(@• 4 ) 





% + f Y » O 

♦ Wy « • 4 ir^ cs*s) 

* 0 

Tfaese ar® siifillar to ttie ¥»3 l 1 Imowi botjudary layer eij^tiena 

the length scales in the bomdary layer are different 
and corresponding liwlt is also different# 

there is also another difference between a classical 
boundary layer problem and the present problem# In the 
forwer, there are two eeiocity conditions t® be satisfied 
at ^e solid wall and one at the outer edge of the layer* 

Here there are two velocity conditions at the two walls* ' 
thus there are three boundary conditions in the classical 
boundary layer pioblen and four in the present case* the 
difference cen be accounted readily by noting that the 
pressure gradient is kno^ in ths boundary layer* Sere it 
is not known* If mm eliminates pressure » the order of 
the differential equation increases by one thus accounting 
for additional boundary c^ditlon* 

tfilh the streanfunction giwen by 
f '^x 


Ci«#l 



C2*7) 


•ipftticm C2*i) 


2 n Ss. 

a.8fc 
2«8e 

t»st 

2.4 omii« Of' nmwoB m cmmmKfiQU of soowiois i 

lov lit pmm@d to OQCkstruet tjoalftleal soXmtio&« 
valid in diffomit partt of tti® flow* TIio ioeltioo to ti#« 
aoalftioal tolotioos m» ‘based on sovoval oontidoiratlons* 

It was antielpatad tfeat as tbe initial pfofil# is not 
analftio f soma singnlavitios will fasult and tlia nsnai 
nmarioal sotioifias wotiid not be foliabla* 40isO| sSnoe mamX 
tnaaerieal •ebtemas are mt designed fot tlie revezwtd flowt 
sone naiov i^difleations would be ne^iired* . flit fesnlts 
in analftioal fom wotad be valid fov mt atwbile nnwifleal 
feault tfOilLd be valid f&w speeifio valnes of a* 

file outline of tbo «ietbod is to eonstmet solmtiffiai 
for Imvpt 1 and alM for seall X and to nateb on pateli at a 


wheire Pjf CX> is the pressure gradient wtiich depends on 
X only* The boundary conditions m t%ms of ^ are 

X * 0 I V <S/2) [Cf/a5 - © s<(¥/a/^ 1 

*» 1 0 <a$)lj^l 

T«tli?'*tl f » 0 © ^ X < CO 

lt«0 |kwO X *0 0 ^ X < CO 

^ tt 

The flow approaches a fully developed velocity profile 



suitab3.« point* 


fli« «oa«truetlo» of solution for i»»li X turns 
out to to fO*T oo»pi®K a®<3 now ideas had to to d®v® loped 
to i«e®t the ocwploxity* Thro® different series solutions 
eaeh rmlid la ©ore, or layer, or tuttl® are roi|uired* 

They has® to to iiatehod m that they are wutually oonpatitle 
Although watching of parafieter expansions is uelllmoun this 
is protahly the first esswple i*i®re wstchJing of coordinate 
expansions is respired* fhis is taka® up in Chapter S* 

For lar^ X the flow way he represented hy an 
asywptotic expansion idiose gauge fimctions depend on X, 
such that as X«^<o the flow attains the fully developed 
■velocity iistrlhution* The treatment for large X has been 
done in Chapter 4* 

Investigation of early stages of the shear laier 
led to a new f awiiy of viscous layers idilch way ha tersied 
as weak wlxlng layers* Weak wixing layers are viscous , 
layers across idiich sharp changes in vortlelty takes place ' 
hut not in velocity* As the suhproblew of weak wiring 
layers provides some information about the mpmsimm 
to be meed later, it is tAk«n up first 'la Chapter ® whieli 
is couple te in itself* 



3. mm KUciBG 


3.1 IHTBODtJCTIOli I 

Tlilii laytr# aerosg wtiich larg» ctiangas 'In '»»loa4%y 
mmT and t^ieh ate distant fro« i©lld walls liaw tssn 
sxtsfisiwlr stadlsd. Caossnhsad 1363 , Loek i@Si, Qi^rtier 
1942) » Mow9'v%Tf s©**wriiat siwllay laysi’s acsoss witicli larg# 
tak« placs Hi worticlty tout not in ■ irsiocity d© not 
smm to tiaiw hem inwstigatad. Hixlng layitrs of th®#© two 
ty|5®s i««y he ©allodl ®s stfong and awak i^ilxing layors* 

Wiiil® m&fmetm diffusion is th® ^tain featuf® of 
th® stnong mixing layors, diffusion of womoity Is the 
priwaiT foatair® of th# vaak mixing lay»ir. Th® woohanis* 
of diffusion may h« th# motions at iw>l»«uiar seal® whifidi 
af« d®so3flh®d phAontMi^alogioally hy yisoosity ox it may h® 
laxg® soal® «ddying motions foiand in tmitniisstt flows* 

Th® oontiMits of this ehaptex thiw lii^ht an th® 
tyy® of oxtansSons to b« ns«d in ih« mixing layor of ^s 
flow thfongh snddon oxpaBsion* Som® information whioh is 
wt portlnisit to th® fiurv through sitddtn axpansion is also 
aidod fox th® Side® of ooiMgilotonoss* 



15« by a suitabi# eddy tlscoaity* 

analysis can b® carried oat aasantialiy along tit® i«s# 
llnai at the anaXysia in tb® Xam^lnar cas«« 

s.a «« MiiiHs Mmn i 


Conaidar tiia staaiy ttfo iiaafiaionai flow of a 
lewtouiaa flaid of comatant proportlas* fba homM&rf 
Xayar approximations are asammad to bo mild and prossuro 
Is takon to bo mlfon as in tbo easo of stfong mixing 
layor* Tho flov is tbon goiromod by tbo o<$natio&s« 


m ♦ O * ® 

it y 

♦ iftty « i) n^ 

and it is siibjoetod to tlio conditions 




* • 0 I 

» > ® I 


y 0 

ti * jcx^y 


y <^ © 

m • Xl^y 

C®*Slb} 

♦ CO 

iy«0». -Q.^ 


y«^ • CO 

-0,2. 

C3«td) 


ioso X and y Cartoslan co«(ovdinatos In atroamijio and 


XI, itwpowid wrtS«ity 

S’?!*. Tli«i» is m tt»«©iitiamity in 
% iaitlaJL saetlon CXLilii^a }* 
flow at tiia iaitial factior*, X^i^ 0 


s«al«» lik« faiplwiE that 2 p ♦ i| 

Ji OB#* Coiifiitlima (S*2®) ana C3»Si| furtbar ro^air# l^at 

f i® #<pftl to i|* fbiis til# ioBstb mSi iraioeit^r temlos t#li«## 
. 1/S 

lilE# m • 


lot til# atv### fonetioB ^ to iiwwii tf 

. 2/3 m 

Y' mjri ta’i'w^Xl ) ft ri). 

1 1 I 


— ■•l/i 

^ • It X) */i ja,^l : y 




fUm tlio •^tiofk of mtibm toooifioo 

f ♦?? - 1 / 8 ? -0 


vlWM srliM dmetta dtffaMntlatlaB* Th* eontittona (t*ae) 


Tile third hoimdary condition required for the third order 
equation (3 ,4a) locates iMm layer. These are sereral 
options for the third boundary condition. Two simple 
alternatives are 

? ( 0 ) 0 

?" (0) * 1/2 {l+Xl/Xi 

If f is a solution of equation (3.4a) then 
af (a’^+ h) also satisfies equation (3,4a) for any constants 
a and h# Since one solution *f (^ ) satisfying equation 
(3,4) can generate a family of solutions f ('^ + h)| the 
third boundary condition selects one member of the family* 

Also the above transformation is of some use in direct 
numerical computation. 

The weak mixing layers having a similar profile 
are essentially governed by one parameter describing the 
ratio of vorticitles on the edges of the layer | as shown 
by equation (3,4)^ The parameter enters in the boundary 
condition but not in the equation. It is possible to 
reformulate the problem so that ,the parameter occurs in 
the equation and not in the bomdary condition* let 

r 1/S . t 

f ('^ ) • [s -a^/ ijCL^*JCl^} J f(‘^)*l/2 IX CS,6a) 

_ _ -I/s 

^ «» [s IZ J ■ ^ (S, 6 b) 

^ « (.a *n. mx! -XI ) 

1 8 I S 


(3.5a) 

CS,5b) 


(S.dc) 



til# #^atio£is (3*4iL^y C3!»4b)y C3«$a) and C3«Sl3> %m&m» 
f*** * if * 1/2 X>]^.) A I - 1/2 (f* 4 ^%^ f « 0 (3*7a> 
^ 2 • 2f" jfe 1 

f m m 0 (3*8a) 

f*(0) ** © . I , ■ 

fii« pafasietw A/i is the ratio of swrago torticity 
of tli® two stroavs stod tho wrticity change across the 
mixing layer* f** ) ]Oipre»ent» to stretched ce-oMtoates 

the iepartnre of local eortlcity from Its average value* 

: In liie ahove formulation _Q. "wm asswied to be ■ 

1 

greater than -Q. ^ * itoemlse a sll^t modification is 

2 

re<|uired in C3*#)| bmt toe problem is essentially the SMe^ 
as given by C3*7)* 

in analytical itetood of solution of (3*7) is to find 
solutions valid for large positive f large negative 
moA .sstall and patch up toese solutimui' ' at suitable 
internfediate points# these solutions are 
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' i 






f®iP' 


C. 


m m ^ ^ r® C ll 

♦ (Sit#!) (©a+S) \^2] ^ fijy «>6* fo? «< 0 

• - 4 ( ^ -■ 1) [8 {■ ♦ 1) ( A «lJ ®* - 1 

• O* * 1) (Sti ♦ 8) B_ -7 I for « > 0 , W.lfflb) 

W • SJ 

B * 0 f#r « < ® 

« 

» j^ (a ♦ Cii ♦ i« ^ ^ * L * 


im « I> CiB^I 


*1 1 

hao ‘- 


1/2 in ♦ iM» • a ♦ U* i»*10el 


Vl Wl - ♦ «> Vs 


If tlia »a|iiitiiS« ©r wrtlelty ©f tfe© ijw«r 

if iawf mm mm ©f lawt# ftiafra lund if rnmrnm 

: . . . ■ ■ 

naif if A liff ia mm tmm CSfll# ^ •^f© 

•©latifa if fpplitfteif ia iliif laagf itofa A« i« 

fliif ffif lortfff i«iffi « ©©ftiffl rnmm mmn la f©»tf®t 

fita »%»ti«sffr nais* fa®© cs*®ai •»« <i*iss> mi© uniJMifi % 

f r- ^ \ ^ ^ t 



If 


hi-} 

B « - r®0 5 Va 

® L- 

fayw>l C3.10«l) 

fh« recursive relations express all the coeff icieiits 
All., «ft4 C., to terms of lower coefficients a»d to all 

M W * 

seven coefficl^tSi nmmXj^ A^ i • % • \ t t f 

Mia S ( or to esse A - 1) are left to be apeelfi«J. 

in® of these coefficients is determined from the aititrarlly 
chosen third homdary condition (S«Sa) or (S«Sh)« fhe 
rematotog six coefficients can he determined hy patching 
at two sultahle points. One patching condition is that 
f and its derivatives are eonttouous* An alternate 
condition would he to minimise a positive definite measure 
of the Jumps to f and its first two derivatives. 

1«B tOlBOLrat IfEAK MIXIHQ LAXEI t 

Consider statiosuiry tui^mlent flow. / It is assisted 
toat if the l«agtli end velocity scales of mean mcticn 
hehave lllEe wf and effective Itoematie viscosity 

hehaves IUk# § as a first appmimation* flie eons^ 

tsnt % is independent of howdary ff»d toAtta i oenditiixns* 
fhe ecpations of mm moti^ mM continuity and houndavy 
condittims are similar to and fhe todCK p 

ia found to he m» ttm eQuatton fltlK As In toe laMinar 
caset the conditlmi Cf.lc) and Cf«fdJ re^rt that y and t 
are e^ial* lat toe meisi streaii fwcttony^he giiMi hy 

Y'eiahi®? t ' li.llal 



mik %»% 


m 

m ' 

fii® ©fsattoB of mm ^tioo reduces to {3#4a| and feowadary 
©ondttioBS to C3»41t5), If li sad 1 ara takoa to bo K/2 aiid 
-Q.jK/ 2* Thus th9 similar prof 11® of turbulont wak 
ffllxiag layer i» govemod by tb® sa^^e ©QUatloa a* it# 
la^iaar eomtarpart* : 

S*4 BSOTITS mu DlSCtJSSIOi t 

fli« result# gt-wm la fable 1 eere obtained by direct 
aiiserical integration of C3*4aJ ifitfe boundary condition# , 
C3«4b) idid C3«Sa) and bare been sheen in figures 1 to 3* 

This three point bomdary ralne problem ©rer doubly infinite 
range liras solved by iterative schevte* f andf ' CO) 

vere assitfned and integration was carried omt for positive 
«.a 5) . Ih. InltUl «1«. ?’(0) f (0) «»>• 

adittsted till the bomdary conditions (B»4b) mm satisfied* 
Calcniations ’were pevfomed on IBf7044 of the Indian Institute 
of feelwologyt Kanpur maing fonry^ ordar Banga^Khtta *9elliod 
eith dills variation* Steps of ^ • 0*01 aare found W be 

*.b.ek«lth..tepofA7.0.0Mf** 

Shoaed that the fonrth digit 1® vortiolty 
Has' not affected by the step else* 

' Idle gmmWk of aesh and itnmg nlxini layoirs is ^hite 
diffS'fatit alMm flanr is iHsinar* the tiiiclsiiess 'Of tte' 



tl 


laytrt liic« x*' and x 'Mmmmw 

th® grcurtli ©f lay««5 of tot^ t;^® ii llBoaap idioji tli« fiow 
is 

Figiar® 1 shows the velocity profile which nm fee 
deserlfeed m a roaaded comer profll®* .Figure 2 shoiis 
the varlatioa of vorticity (to the feouadary layer approxi* 
*natio»K In the central portion the curve Is fairly atraight* 
file change of slope is quite rapid at the upper edge and 
relatively less mpld at i^e lower edge* Jyi the vortioity 

diminishes the profile developes a long 
tail towards the lotrer edge and the thiohness of weak 
mixing layer increases* 

Figure 3 shows the distrifeution of normal velocity* 

It has a linearly inereasinf feehaviour at the upper and lower 

. |i 0* this fe^avlott? 

1 

is different from that in strong nixing layer* fo see how 
such a feehavlonr arises | suppose that the nonal velocity 
approaches a li«it for large positive and also for lai^ 
negative ^ or increases less rapidly than ^ * dny •treanliM 
suff iciwBtly far away fron the mixing layer then has aero slope 
as v/m approaches aero* Row considor a control volisso hounded 
fey one susii stroaniine at the upper edge and mm at Ihe lower 
edge I and fey two sections x « and x » a^* Since ^e flew 
accelerates in the downstream directlMf and ainoe the areas 
at the upstreaa! and dewnstresm sactions are equal t the xaa# 
efflux throu#! the downstream section will CKoeed the mass 


edge of tfeo «lxin.g layer idien -cx^jCl 




iiiflm til® i®etioo. Tiiif ’riolatioii of 

eofitinuity sliowis that th« lioiwal wloeity caimot &ppm&da 

« Iteltlng »alu* M ^ 00 Mia - CO. Ih. .fco« 

airgixiftfmt eanmt ho appilad to eas® iihain th® fluid os 
th* lomv ®%« ia stationary (Xig/jn^« 0)* 

If hko third touadary coaditlos C3*S5 i® chosiKi Is 
a diff®r®nt wayi th® shap® of th® u valoelty prof il« 
remains unaltered and it merely gets displaced* Boweiery 
the shape of the normal yeloelty profile ehangesi as the 
BOiwal velocity changes froai ^-^f *to f*t>f under the 
tr«aafoi«atios f )■ to f Cosseipwatly, for on® 

particular Ohoioe of housdary condition one can ohtais 

a normal velocity pprofile is nhicdi ▼ approaches a llmitisg 
value at the upper or the lover edge* the previous argUiMest 
however shows that the normal velocity cannot he made to 
approach a limit for ^ ^ oo and at the same time for ^ •*^-00 

The weah miacing layers for idiieh ^ i# positive 

are <iualitatively different from those studied here* The 
Initial profile in this case has reversed flow and hence 
the ccdveotivi influences are carried In the negative 
s*direoti^* This upstremt propagation would also rsnder 
the ipecificatiiui of initial ctModition with a sharp comer 
unrealistic* 

a icneralisatiim of the Initial condition C2*ta) end 

would he 


s « ITq+j::? <S*lSa) 

mwiQ-^xzgy C3.i2b) 

Til® pwil>l«i?i is !ioif«*9»!f to»i«iiiiat l0 i» mt 

s#3p®* I'S 1^1® pifft'vj.ottf Pafi® ! mm 0 I y®X0Cl.^y la .ite® 

Tiseoa# i®y»r was o(i) aai la tli® sscoai ®ai« %% 1® 

If €»a® tak«s 

( 3 * 33 ) 

th« goir«»iag «Qiiatloa l^seonss 

^0 ^**1* * ^ t3*1.4) 

whsa® isli® tiasaswsas® JL®ag%ii s®sJL® ^ Is ttslosa %© fc® 


4. FOB umm x 


It if. asfiwfd that at downstream* far fro?" the sudd«a 
expaBSioa th® weiocity distrihution tends to fully developed, 
flow.* Th* taEpauslosi for straa«»’ funotion for lari* X ostu 
h« written, at 




t4*l> 


idi«r® ^ as X OD and * 1* /q representi th* 

str*a«» ftmotlou mrvmpm&lMi^ to fully d®y®lop*d flow* 


*» i 3T-1® >/a 


Substitution of (4*1) into (E»?) yields 

^ »,/ 


(4.2) 


2 [<^ ♦ 2 «>. < '!r >1“ 'k 

” If w to.t g* g, - -x„ g„,, g, 

a# f ** 1. Hone* g • «sp C* X X)* Also • A g 

0 ^ 'Vo- r*s Vs y+s 

This oan h* satisfi«d by tahing g^ « «aqp> (-r)vX)* Mith this 
ohoic* of gjp I a<ld*tioii (4*S) reduces to 


for all »• FjjCX) is atstwed to be 2*^, % uh»*« 

are unspeoifiei* Ilewlnatini % s frow (4*4) leads to 


C4*4) 



S>6 

fSi© t'omdary eoMitloni on 1^ i«pl^ tliat 

f • t t I ■ jit « » © t » # © 

for tHo oas 0 ^on m is e<|^l. to me (4»5a9 l>) is 

ji'V xp3/i;ci-'y^)<'+3)^j .0 

I • * 1 » « y; » 0 

«ti«r« X appoars as aigaaialua* 


|4»Sa) 


C4*SI>J 


■(4.6ai 

C4*6li| 


A siitplo Mtliod of dotoniimiag A is to sooIe: a soriot 
soiutioa of 1^0 fof« 

; ,# 

^ • Ik “^a <*•’«} 

idiorci thtt oooffioionts ^hdopoad oa X « ttia rotursifs 
roiatioa for ttio oooff ioioats is ottainoi frost o^ctiatiim 
C4*6a) as 

Ca4«) «a^) o< ^ 


cK « « 9 j^^X 

(ii 44> la#®) Ca#®) 


C4.Tfe) 


for all &• 


4 ioaorai soXittioii of C4»€a) «aa 'Ni mrittaa as 

V-, - ^ '^’11 + 'V'15. /^-i Via + P<i V'n. (4.8 

«I>*M 54 . '■ an ItBMUPlr lna*pM)*«n» funatlao* nf * 

tit*: „ ■ 



fatisf^ing C4*6&) ant s sto 6on.stan.ts to to oiralmato4« 


file function# ar® ottaiaed fnana ofaation C4*?a, t) 
on taking ^ 0 * » if » 0 wi gat 

X/m ♦ )? /1JJ20 /- C4.iaJ 

taking » ^2 • Of » i 

* Z X/40 / -► A/140 C l^SX /8 ) X*^-44,0l>) 

taking ^o* if • ^2 a ^^3 * 0 , 

» 1 * X ^ / • •*• .14,9®) 

ang taking • it » 0| • if ® 

)^,4« i • C4*M) 

Will® obtaining tha atove solution# tka coafficiwati 
^6 I ^1 f ^ 2 . ang ^3 an® ©feosan in a way sack tfeat tka 
oW an4 awan aoiution# i© not gat ^liitaa* 'Vn aft4 %j_ mm 
o4d functions of X wbairaaa yf^ aii4 V-h afa atan* 

For a noatriwiai solution eongltion C4,01») wouli 
rofiidfa tliat ‘ '■' — ' 





%(■>> 




hU'} 

V-iitrl) 


r,j(-o 


^ilC-o 

V, U-i) 


Y,^(-0 



uliieli oia siiflpiiflcatlon yiatlSs 


^ hiM 1 1 5 

L 'f'lihO Hil+'/j L 



®if ^ CX) » < X 3 X % ( X ) « 0 


( 4 . 10 ) 

(4.11) 


^ %m first m& mmnA mrlf Iraeksts ia 
si|ii»tioa (4*10) r®»p®etlir@ly# 


fli® s@ros of Q(X) are the elgim¥el«» of (4*6). 

Q( X ) la (4*11) has beea rep roseate! a« pw^aot of <i^( X ) 

). n» *.«.• Of Qj^C X ) M. th. 
eorrespoaSlai to o4d olg^foaetloa aad that of X ) to 
Ofoa eiioiafiaetiiia* & almple oo*npMtatioa fiolde ^e 
eigeairalao. (fable @1). 

dillie aad Braadt (1064) atteitpted a direet aiiaorioal 
solatioa of the eqjuatioa obtalaed fron the toe(plete latter* 
Stokee e^atioa* ttiejr are ooacemed id.tli detelopiiieat of 
Foiioiiille flow* The reletaat resalte are eoataiaed ia 
Table SB* aillie and Brandt are ecmeefned oaljr tdlh odd 
eigenftaietioa and find that there is a eritieal leTnold 
niMber hbott idii^ ^ eigimtalno is leal* (B^^. « t*4d 

•Wlp* 

Wilson (lfi») again starting ftoe eewplete latter^ 
itoiees e^ttaiion finds that for large t there are two 
•efoenaes of eigentmlnes* Both se^foenees are aifttptotioaUir 



msX ai B m * Tlit wawtoar of thm first loiponc® 
ar# 0(E) at H *-4- cd* . fiia g®eoHd sefoanc® eoniists of: 
oigon-falmot art 0(1) ms 1 oo* tha flow win t>e 

ioainatai Tsy the aaatar of tiio soooad se^ponc® for largo 
1 and tho disturbane* will porslat downatroa*?! a distmnc# 
Oil)* Wilton uted a ctpansion in powrs of for tbo 
oigonraliio* 

^ « In ( (4*12) 

fb® r®sultt ara glwon in fabl® SO* CosipaPiton of tlio 
prosont rosnlts wiHi thoso of Wilson and OlllSs and Brandt 
sbows tbat tbo Unit oonsidered in Cliaptor t md the 
oacpansion (4*1) ut«d h«r® prorldo eaooliont approxinationa* 

Intorostinglf Hi® dominant nod® of disHizbane® 
wiooity if antlsywnotrloal- It follows at tMs will 
liaro fowor »i>os tban ftymotrioftl woloeity poniii1>ation« 
Booaua® of tli® atsnatod ayanotir of flov In tb® problM 
of tnddoa aotpanaion tli® oigonwaliios of our intoroit aio 
odd* Tborofor® «» got 

^ m • 0 

h - A ( ^11 ^ 

!*•»• -r. - 

flio liagnitudo of /6^ bat to b« d«t«r«lii«d bjr patfiMng* 

fbo problom of finding bi|^®r K. It ^ain^l®* for 
li>/i tb« bitid Sid® of C4,#a) it a forolng fwotion 

Wbleb it fttowii* 



aisiswptioR ®j| • ia p0i®ll;Jy® to gmrnmllm 

ttit* In %fe« following 


K" • )t> ♦ 2 oac]^ C • n A X) % If) C 4 

wt^ow tlio 4ooll# ®«rl#« is soitally orrang®^ *o that tho 


giiitgo fianotions «ro stieootsiTOlj of htguor ortor, and 
Aid < Amf/ ♦ A^ or® tllto oifoowloo of 


y-Z + ^ m [(i/iX I- Y V y-z -I- 3 Kr,,]^ 0 {4-16*) 

r • ♦ i 1 )J;,''* 0 ( 4 .uit} 


f&o oigMWaloo® stoptisis f ton ono* fli« pfobloft 

of flMtof t , y\>i i* ilffiilaip to lisot hm loon dtoowsod 

fyfST) 

ftorlior* 


' tno oifiminiloit ifi^ nofotivo tool i>«rt «i«o 
i»»o«. € Wlloon IW >* fo» low » tn® only oigoBomln®* 
iilio*® tool fort l» nogatiw «»• ^o»» OCa)* It oan tnortfor® 
to taii tl&ot tto diotnrtoneo la^stroon of tlio stiiSon ogfiiiaioii 
fonotrotoo • diotwoo lain®® ta® rolowant Bofiiol^ 

MBiitor fot Of itfoow ^Hbmoml i« IUbo tiio ftop^ia noiilor for 

lUPwit 




mLPm&im fob small x 


IStiOSTOlOi I, 

In tliiB eti»i5t«r solutions for various w^gions at outlinad 
la saotloa 2*4 &m oonitfuctai and watchlaf is contiioistd in 
th« llwit X 0. A gouaral ■"tethodolofy to construot th& 
sonios solution nould T&o to Krrlt© th® indspandont ifmilatiat 
jilt ttia F o 3!^t' 

5=^xr. 

nhan A M Sv,?:*’ * Hio iraapiatlat 

aro to ©hootan to tak® ©ar® of fractional pounrs of X by 
taking It appr©prlat®ly* fk® dot laitlon of ttotns that X 
it atr®t©h«d about th® ©mrvo A(^) t. strotohing funotion 

b«in* 

' tha form of mtmm fmotion it ©iiooton at ' 

^ 2ki C ^ *+■ CS*lb) 

It it axp®ottd tliat log t«f!fi uHl also appoar in ten® imt$m 
at sofRoidiat lator stag® and to a gonaral fora of stvoan 
fuaotion it vrittiBi. shall bo t«k«i taro imlati ra^piiapod* 

By tpooifying kf tho loading ton in t!i# o*p«»sion of 
atroaii funotloa ©aa to adjuttod, sabstitution of to 
12*7) loads to 



Cs.iii) 


k+Y + ^^r(^+m 

+ 9 c. Jp ~ Gsv^s l-ci nt. 

Jf! ^ ^ ^ 

% +2^<l 3^[L<+r^^t g J + ^/3r 

k+T+& + ^« M+vn 

-9^[(k+s-i4'+9/]-Ck+s-<)^^} -Gsyv.W. <«- 2 t) 

3 r 9^' -(k4.s-i!j 9;g ' f G^-f«v.k js. 8,3 

fh«8® eon«tar)ts art related to preasure gradlmt* 

S ,2 SOlOTlOi FOB IIE COBB BSQIOi I 

la til® oor» m tak» m » 1. Slao® tb® »tr»«st fmett®® 
is to matoti til® supplied initial oonditioa at X » Oi l ia 
•pooifiad *®fO* Str®aa funetioa i« mdt%m m 

f-= <»•») 

vlWM 1«T/«<X), «(X)» Th\,3^* *10“ *• 

On au'toatltutioB (s*t) r«ilu®«a to 

T 4 -S -4 Xtfc 

i^«r® i* oonatant flwa "bf 

m 


k+V 4 -£+^^^ h-i-vn 



m 


til® |©»«t is for » ♦ 1 » 0 

T+sf i "=0 

ai>4 i» satisfied for «r?y F^* ^h» iulti®! eoif»diti©ii 

re^Mirts tin at 


• C S f • r ) / ^ 

0 


so that %* • t 3/18 M i - f } 


{$*m} 


k fpnarml fo» of th« aciiiaticm (S*4| goraming F, 

I 

©an eoiiv«ni«»tl.F ^ written dam ^ 

m ( Ro^ - Po^ Ri) = n C Pb ) _ C I)^, ~ <S< 

4 - ^ tt ( Fr FK\~r Fn-r} -t 2 S-t ^ TT /v FJ- 




j£»l 

T4S+i*h 




irhere 1 L«IL/%» 1 ^/« 

^ ' li li ' H IBI 


Squatioil <S#7) is eiispilar as F^i tte vtmtfMmt of 
fwlshes at f • i 1* its® teeaas® »f tlie assuwed sywwetty 
the eowpiifierstaiy somtiiiii wlii he »e*e setd ®«ily j^ttieiilar 
Aateitai wlli watilhiite t© the soiatiea. 

f®r tt • I einatseii <i#f) ie 

( i|^"i^ > - f - < I}, . » / . w^> 

Its mlmtim is ' . ' ^ ^ 


f , • tJ 


i/t) U-l®> tl • MS,/*) >*• 

«• ,:■'♦■ " ■ . ■ ■ ^ ■ IfIS 


(S.Sk} 



tl'l' 

wtitr® ^ i f 0 ^0 ) / a i# coBstmiit# 

r'^ « (3%/25 it^t ) ^ %Qg Ci'^i/i-i) ♦ t| (s*se) 

fli® abow tom of f indicates a log singularity at Y » ♦ 1* 

This sittgtilarity can fee attributed to the fact that at the 
edges of the cote the condition that Telocity is aero (before 
the sudden eotpanslcm) has be« relaxed and con*e)|a®fitly the 
local behaviour is different than that of elseiid»r«' In the 
core, fiance a different kind of expansion is needed to 
explore this regloiit nmely the shear layer* It is anticipated 
that '^e sin|ularity uhich shall arise ifill be mild in nature* 

mmok FOE THE @11141 lAlEE I 

The shear layer can be taken as a special case of 
weak wixing layer idiere a vortical #tre«i cowes in contact 
with fluid at rest. Baaed on earlier InvestigatS^Mia (Chapter 3) 
we choose In equation C®*S.J § « • S and k » i and write the ! 

independent variahles 1 

[ Y- (o.#,) I 

i 

Where ^ * iCXl and the sthetcMng fmcti<« hC Is 

■■ ■ 'i 

simply^ h|l • the fore ®f the strea» function is ehoos^i as 

> = 2^ f [I05 5 . SJV^ + ^v] I 

$1 ■ '■ . ■ ■' ■ ■. . ■ , ' ' , ' ^ 

Iquattcms CB*8> reduce t® . 

=: ecmsfawr 


(itlCa) 


»'<. Z [0i f(r+^.;^^+94+ 

V4S- « 

-' ^(e>i/}Ql-a^f <§#1.0fe) 

I^ol [( Tt^) sX- s: 3 i J-' , 

T!io left hand ild® of (S#10) shouid ht indopandtant of ^ md 
thaao oouftiffAtB aro roiated to p3?««inr® gradient* 

Though a ganoral fow of the governing eipatioaa 1« 

written the proceduie i» that g la t»km to ho aero rnilosa 

11 

need arises* 

In the lowest order e^nationi the ohoiee of m has 
led to the ap|>earenee of Inertia and diffiutiim terms together* 
fhe lowest order efoation is C is tahen sero ) 

K( ^ ^ KK > ^ ewsstant (S*n) 

this equation has he^ well investigated In Chapter 3* 

for large t difftisiett term drops t© sero faster than the 

eonveetive term* Also f h^iaves as a foadratio In pXm m 

exponentiallr small term < SSI ) for large 1'^/ towards the 

vortieal stream and as a eonstant pins S8® towards the fluid 

h 

at rest* this ®an readily he seen as an •nbatitntien 
in Ci*ll> one ohtalna » • « or ©♦ Writing 

~ « ,if+ ^x-Tj + QST , ^ IflTge /’] 1 W.M) 



IjIi® hand aid# of C6#H) is i^C4a^i^ •• a^) nhioh is oonftani 
aiid hmm th%m is no rsstiPiction ^ $ »2 *s* 


Th« noxt order gOTreming elation Is ( Is taken a^ro ) 


+ li^( Sf”f ^ ♦ 2f^f” 

jx * 1 or 3* Let 


faking f <v 

w 


1 


*• * eonstant ■ CS,13) 

I , 

^ for large n\ m find that 


for large I’] 1 * For <§*X35 to be todapendent of -n are gat 

' E ^ 

ettliar ^ «g • ^ ® or a^ * a^ * 0* M« assisna the first 

one* If it is not true then left hand side of (5«ii) is 
nonsero and henoe the leading tern in idie mpansion of pressure 
gradiisttt will be 0(5 * )• get therefore fron (5*13) ■ 

• C ®V®*1 ^ H CS.lSa) 

with additional oondition 

6a^ ^0 ^ ^*^*6 **" ®*x*V ^ 

Httlle earrfing out a aiailar analysis for the wUsiai layer 
for large '>] towards tJae fluid at rest we obtain 

fo C®*lda) 

b^.^] t b4 i-eST 


Ci*l€b) 



th® left hand ©Id# of Cd*ll) and (©•IS) shall he »®ro 
and with no additional condition on constants hj^ *• 

Higher order equations can he analysed in a si»>tSlar 
Twafif-er* shall return to it idionewr th® need ariaes* 

S.4 irn’CHild OF COBE aHO SHlAE h&fES EXPASSIOISS I 

A general watching procedure C Kaplun 196? ) is outlined 

helow* 

The irarlahles In core are X, f » J/MiX} and the streaw 

1 y) 

function eicpansion Y- ^ Z,^ X F^C?)# The expansion for velocity 
valid in core region is 

S m • y f * / mm 

n 

■ r’ ♦ X ( r’ - If' ) ♦ o«®) (b.W) 

W 1 1 © 

1/8 

In shear layer the imrlahles are q » ISKJ and 
^ The streaw fimction expansion Is 

- 5n ^'-n J * Therefore the expansion 

for veloelty valid in shear layer is 

2.f n , / n / 

9«+'f^JAo^ , 

Introduce the intermediate variahle 

Y|u = [\~ H()()]/jUCx) «•»> 

mwM «x)3A^ttO i»«). nut u, 

is m intemediate length scale Cit does not repieswt 
ooeffloiwftt of visooiityl* Therefore 


’]= Y|U (M/ko£) 

V ' + Yju (A</tt) 


Tlmii tli« mtclilRg coMitiOE it writtta at 

- cS(oj] ! =0 («.») 

la th% a’bo'w li’Rlt, X -VO far fixed f,\>* Th® watclitag 
it said to bo doa® to th® ord«r of tit® gauge ftmetloa 



Coatidor tli® it»at«hinf 

'^’r^ [G(u)-g(.o ;]/4 ='^ ( 6 .si) 

FdY'?)= Ci/a)ti-Y'^) 

= (c-V2j[2Y|j M/h+V C^'/w] 

J-t) C'>)) - 

~ 2.«iY)ui bro-e '*) 

E 

Ih* nktehlni laad* to mj • *• and ag » 0. 

Ihaftfot. ig - 0 , .g - 0 . »«w «- a.trt»te* 

/!■«) [S(Uj-<?(U)l /£^=0 

(M /ko£j + Y 

E E ' 

•ad satiMBf i««df to % • • a*^ ^ 



C5,23) 


of gjg* tlio ©quatioii go’fomiag it 

fli«r«fO'» i <v‘^ * Lot 
2 I 

hTj-g-e />^) (S.84) 

and tlion loft liaiwl aide ®f C®*83) 1» ll3co 8li^«-a. Ii«ne# 

0 2 1 

thoro it no rtttxlotioii on tli« oonitmts i»| 't» ftm tlit 
liatoliiiig 

fcttJ) ~ SCu;]/ f 'a o CS.2§> 

n 

m g«t tliat if log log^ *4^ ©t fhim 

r«tt»letioii o» jif mm» that thar® It a narmw ©wrlap 
(••g* fu ^ ^ )• Othatwif® tha niatahlni e«a ha 

acecwpXl Aai only liian and a*» aato# 


« « Hftfer, f ® A 

«9 • 



m V go«i to »©ro mpmimtl&llf in the ilioar lajror* It 
can therefor® t® thought that there is a ImjitT cios® to the 
wall X * called wall lay®r$ idiich supplies the entraimwat 
corresponding to the leading tern* fhe wall is considered 
a singular surface as our assumption (that' the length scale 
in the direction of channel hehares liJce K) definitely will 
not be valid there# Any ai^itrary surface i%o the order 1) 
will be reduced to a plane surface in the limit R -v oo# 

Hopefully this difficulty can be bypassed by matching 
Me therefore retain the viscous term also* 

r- ♦ a-.) ( f . f ;r; ) - 


By redefining the i»iep<«adent variable wad tiie #®ostant 
w® write the Above eijuatlou as . 



m 


t!s@ tiOMidary coMitlou* FqCO) « F^CO) » 0, * 

Atswiiig a ierte# solatioa 


r-> ^ 

-- Iv, n.v 


th@ racurslw ralatloa for th® coefflci®ntf’|^ar« 


r, 


V>-^4 


[(jof 4) (Ti+'VI-Jn't V Oo-^O] 


whem r'o - Hf 2/ Ao /4 # flw aatooi®. 

i» to t# found hy matotiljif vltli t^o shtar 


In thm9 oaotions a Kotliod of oonatruoting solution 
for mall 1 has h%m outlinad and soaa staps haw hoan 
eaafylad out* 



6. COliCi.0SXOI 


A profciskt* ol* 'thd cXitss lu^atiiBal 

separated, flow has iawstlgatai for large Seinoias 
nw»lBm using an aspiptotia i^etfioa. fhe pmhlm ehaasm 
is the flow in a sudden expansion in a channel* It is 
assu'^^ed that the length scales are O^IIJ snd OClJ in stresw* 
laim and transverse direction* this is co^patihle with the 
computation r»ade hy Hung and Macagno (l©66)* later Maeagno 
and Hung (1967) atte’spted numerically and experiiswtally 
flow past a conduit expansion* their findings ^ee again 
support the haeie assumption* 

the main Constance of the assunption is 'that ^e 
iavier^Stohes et^ations tecome in this limit homdasy layer 
equations* the proMem differs In two vi^ ffom a classical 
boundary layer pfohlew* first t there are four boundary 
conditions instead of three* Secondi pressure ilstrlbutitwt 
Is not Icnown and Its 4Nitefwilnatl*Wi Is a part of the ptobleis* 

fho mothod of aolutlon is to constiutt coordinate 
eapansions for small 1 and for largo X* tor the ^^Mislons 
of the first type I one needs three limit e*p<snsio»s ralld 
In core f miring layer and ^le bobble* 

the matehliig procedure has been outlined* Intlthtin^y 
the concept of matohing In problnis of paranetor inrtuibatloa 



flow at largo dlstamaas from the «aSi«. expansion 
(Chapter 4) leads to an eigenralne protlew and the deterwlned 
•Jgeiimlmes are in exeellent eg »e»^*nt with that ^ohtalaed hf 
Gillls and Brandt (1964) and milson Cl96®)| they hare nsed 
the co-*plet« Hairler-Stohes eipation. This that the 

ll^it considered does bring omt the osswsitial features of 


the flow# 


To complete the worlc, one womli need aiditioaal teiws 
in the inner expansioni patehing the expansions at s^^all and 
large and eoatparison with the resnlti of Maca^ and Bung* 
Daf initire Judgmant ahcmt the validity of the expansions m& 
correspondiJii lieiita nan he eade only aftenrafdi;« At the 
present stagef no eajor feature is rerealed whiah would he 
difficult to explain In teiwis of ^nalitatiwe physical arpwie&ta 
and thla is an Indieatiim* hut not a deeonstratKm that the 
essential aspects of argaw^its are calif • 

The study of Initial dewelop*icnt of the mixing layer 
has led to the identification of a family of cisccus layers 
across idiich sharp changes In w^rlicity occur nn^coj^anlcd 
hy sharp changes In rslocity# A complete stisSy of these 
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liiy®r» wmal 5iiter®sting properties • For osca^pl®, 

its thictooas varies like wiiike tfee stroag i*f$xlag 

layer %ih%ch grows like 
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Fig. 2- Vorticity profile 
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Fig.3- Normal velocity profile of weak mixing layer 
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